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Abstract 

The spinless Salpeter equation represents the simplest, and most straightforward, generalization 
of the Schrodinger equation of standard nonrelativistic quantum theory towards the inclusion of 
relativistic kinematics. Moreover, it can be also regarded as a well-defined approximation to the 
Bethe-Salpeter formalism for descriptions of bound states in relativistic quantum field theories. 
The corresponding Hamiltonian is, in contrast to all Schrodinger operators, a nonlocal operator. 
Because of the nonlocality, constructing analytical solutions for such kind of equation of motion 
proves difficult. In view of this, different sophisticated techniques have been developed in order 
to extract rigorous analytical information about these solutions. This review introduces some of 
these methods and compares their significance by application to interactions relevant in physics. 
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1 Bethe-Salpeter Formalism in the 
"Instantaneous Approximation" 

Within quantum field theory, the appropriate framework 
for the description of bound states is the Bethe-Salpeter 
formalism [ 1 ] . Therein, all bound states of two particles 
(in fact, of any two fermionic constituents) are governed 
by the homogeneous Bethe-Salpeter equation. Here we 
are interested in a particular well-defined approximation 
to this formalism, obtained by several simplifying steps: 

1 . The instantaneous approximation, neglecting any 
retardation effect, considers all interactions of the 
(two) bound-state constituents in their static limit. 

2. The additional assumption that all the bound-state 
constituents propagate as free particles with some 
effective mass m yields the Salpeter equation [2]. 

3. A disregard of all of their spin degrees of freedom 
focuses on the treatment of scalar bound particles. 



4. In technical respect, the canonical transformation 
\x (1) 



P 



of position {x) and momentum (p) variables casts 
in the case of particles of equal mass m for a scale 
factor A = 2 this approach into one-particle form. 

(For more details of the derivation, consult, for instance, 
Refs. [3-5] and references therein.) Refraining from the 
nonrelativistic limit, we get the (nonlocal!) Hamiltonian 

H = T + V . (2) 

This operator is composed of the "square-root operator" 
T of the relativistically correct expression for the kinetic 
or free energy of a particle of mass m and momentum p, 



T = T{p) = + w? , (3) 
and a (coordinate-dependent) static interaction potential 

V = V{x) ■ 

frequently, the potential V{x) is assumed to be a central 
potential that depends merely on the radial coordinate r: 

V^V{r), r=\x\. 

The eigenvalue equation for this particular Hamiltonian, 

H\xk)=Ek\xk) . fc = 0,l,2,..., (4) 

defining a complete system of Hilbert-space eigenstates 
\Xk) of H corresponding to its (energy) eigenvalues Ek, 

{Xk\H\xk) 



fc = 0,l,2,. 



is commonly known as the "spinless Salpeter equation.' 



2 The Relativistic Virial Theorem 

Useful general statements about the solutions of explicit 
or implicit eigenvalue equations may be proved with the 
help of virial theorems obtained by generalization [6] of 
the well-known result of nonrelativistic quantum theory. 
(Ref. [7] is a brief review of relativistic virial theorems.) 
For eigenvalue equations of the form @, the derivations 
of such virial theorems can be traced back to the (trivial) 
observation that expectation values taken with respect to 
given eigenstates \xk) of H — or matrix elements taken 
with respect to arbitrary pairs of degenerate eigenstates, 
\Xi) and \xj); of H, i. e., eigenstates satisfying Ei = Ej 
— of the commutators [G, H] of the operator H and any 
other operator G (the domain of which must be assumed 
to contain the domain of H) clearly vanish. Suppressing 
the subscript that enumerates the eigenstates, this means 



{x\[G,H]\x)=Q. 



(5) 



For the symmetrized, self-adjoint generator of dilations, 
G=^{x-p + p-x), (6) 
and H of the form (|2} their commutator [G, H] becomes 



[G,H] 



dT dV ■ 



In this case Eq. (|5} yields the master virial theorem [6,7] 

(7) 



X 



dT 



dv 



X 



this relation expresses the equality of all the expectation 
values of the momentum-space radial derivative of T{p) 
and the (configuration-space) radial derivative of V{x) , 
it produces the specific virial theorem for a particular H. 
For any nonrelativistic (Schrodinger) Hamiltonian, i. e., 

H^Hs=m+^ + V{x) , 
2 m 

Theorem entails, retaining the conventional factor i, 

' X ) = t: ( X 



X 



P 

2 m 



dV 
X ■ — — (x) 



dx 



X 



(8) 



For the semirelativistic "spinless-Salpeter" Hamiltonian 
(|2j, involving the square-root operator of the relativistic 
kinetic energy (|3j, our master virial theorem leads to 



P 



X 



X 



dV 



X 



(9) 



In the nonrelativistic limit m — > oo (i. e., forp^ ^ m^), 
this spinless-Salpeter relativistic virial theorem, Eq. (|9}, 
necessarily reduces to its nonrelativistic counterpart (|8j. 
Similarly, the virial theorem for the Dirac equation [8,9] 
is easily deduced [7] from our master virial theorem 0. 
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3 Bounds to Energy Eigenvalues of 
a Spinless-Salpeter Hamiltonian 

The precise determination of eigenvalues of operators is 
of particular importance for any formulation of quantum 
theory. Unfortunately, for most cases it is not possible to 
determine the point spectrum (the set of all eigenvalues) 
of a given operator analytically. Several powerful tools, 
however, allow to derive analytic bounds to eigenvalues; 
applications of these techniques to the spinless-Salpeter 
operator (13 are reviewed, for instance, inRefs. [10-13]. 

3.1 Minimum-maximum principle 

The theoretical foundation of any derivation of a system 
of rigorous upper bounds to the (isolated) eigenvalues of 
some operator H in Hilbert space and hence the primary 
tool for any localization of the discrete spectrum of H is 
the well-known minimum-maximum principle [14-16]. 
Its precise formulation is based on several prerequisites: 

• Let this operator H be some self-adjoint operator. 

• Assume that this operator is bounded from below. 

• Define the eigenvalues of H, Ek, k — 0,1,2, . . . , 
by the eigenvalue equation, with eigenstates \xk) , 

H\xk) = Eklxk) , fc-0,1,2,.... 

• Let these eigenvalues Ek be ordered, according to 

Eq ^ El ^ E2 ■ ■ ■ ■ 

• Consider only the eigenvalues Ek below the onset 
of the essential spectrum of the above operator H. 

• Restrict all considerations to some d-dimensional 
subspace Dd C 'D{H) of the domain of H. 

Then this theorem asserts that every eigenvalue E^ of H 
— counting multiplicity of degenerate levels — satisfies 



Ek < sup 



for all fc = 0, 1,2, , 



In the case of one-dimensional subspaces, that is, c? = 1, 
the minimum-maximum theorem reduces to Rayleigh 's 
principle: the ground-state eigenvalue Eq of an operator 
H is less than, or equal to, every expectation value of H: 



En < 



Given some operator inequality satisfied by the operator 
H, the minimum-maximum principle may be employed 
to derive, by comparison, upper bounds on the (discrete) 
eigenvalues of H, provided that a few assumptions hold: 



• The operator H, exhibiting all properties required 
by the minimum-maximum principle, is bounded 
from above by some other operator called O, i. e., 
it is subject to an (operator) inequality of the form 

H <0 . 

Applying both the minimum-maximum principle 
and this operator inequality, any eigenvalue Ek of 
H must be bounded from above by the supremum 
of the expectation values of the operator O within 
the (fc+l)-dimensional subspace Dk+i of'D{H): 



Ek 



< 



< 



{Xk\H\xk) 
iXklXk) 

sup , ,| 

sup , ,| 

for all /c = 0,1,2, 



(10) 



• All eigenvalues Ek of O are ordered according to 

Eq ^ El ^ E2 1^ ' ' ' . 

• Every fc-dimensional subspace Dk in the chain of 
inequalities which constitutes Eq. (I10> is spanned 
by the first k eigenvectors of the operator O, or by 
precisely those eigenvectors of O that correspond 
to the first k eigenvalues Eq, Ei, . . . , Ek^i of our 
O. For this case it is very easy to convince oneself 
that the supremum of all expectation values of the 
operator O over the (fc+l)-dimensional subspace 
Dk+i is then identical to the eigenvalue Ek of O: 

WeD,+, mi') 

Consequently, an eigenvalue E'^, fc = 0, 1, 2, . . . , of the 
discrete spectrum of H{< O) is bounded from above by 
the corresponding eigenvalue Ek, k = 0,1,2, ... ,of O: 

Ek < Ek for all k = 0,l,2,... . 

It remains to prove an "appropriate" operator inequality. 
(Summaries of the idea to find bounds by combining the 
minimum-maximum principle with reasonable operator 
inequalities may be found, e. g., in Refs. [12,13, 17,18].) 

3.2 Analytical upper bounds 

3.2.1 The trivial nonrelativistic Schrodinger bound 

The inequality (T — m)^ > expressing nothing but the 
positivity of the square of the operator T—m may be, for 



2 



m > 0, written as an inequality for the kinetic energy T: 



is determined by the effective principal quantum number 



T < m 



2 m 



(The right-hand side is the tangent line to the square root 
in the relativistic kinetic energy T in the point of contact 



N = n, + L + l, n,. = 0, 1,2, . . . . 

Unfortunately, in the Coulomb case the squared bounds 
are above, and thus worse than, the Schrodinger bounds. 



p = 0.) This result proves [17] that H is bounded from 

above by a nonrelativistic Schrodinger Hamiltonian i/g : 3.3 RigorouS Semianalytical upper bound 



2 m 



V 



For a pure Coulomb potential V^(r) = —a/r, the energy 
eigenvalues of the Schrodinger Hamiltonian Hs depend 
only on the principal quantum number n, related to both 
radial and orbital angular-momentum quantum numbers 
by n = nr+^+l, withrir = 0, 1, 2, . . . , £ = 0, 1, 2, . . .: 



3.2.2 A "squared" bound 

A relation between the (semirelativistic) Hamiltonian H 
and a nonrelativistic Schrodinger operator may be found 
[ 1 7] by considering the square of H and by realizing 
that the anticommutator T V+VT of relativistic kinetic 
energy T and potential V generated by the square fulfils 

TV + VT<p^ + V'^ + 2mV, 

as may be shown [17] by inspecting some consequences 
of the positivity of the square of the operator T—m — V: 



< 



+ + TV + VT 



2V^ + 2mV 



With this inequality, the minimum-maximum principle, 
recalled in Subsect. lTTl immediately guarantees that the 
energy eigenvalues Ek of H are bounded from above by 
the square root of the corresponding eigenvalues f q .^ of 
the Schrodinger operator Q, constructed by squaring H: 



Ek<^/S^k: fc = 0,l,2,.... 

For the case of a pure Coulomb potential y(r) — —a/r, 
the operator Q has the same structure as the Schrodinger 
Hamiltonian i/g of Subsect. r3.2.1l with I replaced by an 
effective orbital angular momentum quantum number L 
involving both the usual I and the Coulomb coupling, a: 



L(L + 1) =^(£+l) + a2 , ^ = 0,1,2,. 



(11) 



The set of eigenvalues £q of a "Coulombic" operator Q, 



.JV 



1 



2N^ r 



We regard an energy bound as semianalytical if it can be 
derived by an (in general, numerical) optimization of an 
analytically given expression over a single real variable. 
Taking advantage, as a straightforward generalization of 
the (simple) line of argument sketched in Subsect. r3.2.1l 
of the inequality (T— /i)^ > requiring an arbitrary real 
parameter ^ of mass dimension 1 and obviously holding 
for all self-adjoint T [19] implies, for the kinetic energy. 



T < 



P 



for all ^ > 



This translates [17] to a set of inequalities for H, each of 
these involving a Schrodinger-like Hamiltonian, Hs (m) : 



p 



2fi 



V for all ^ > 



The best "Schrodinger-like" upper bound on any energy 
eigenvalue Ek of H is then provided by the minimum of 
the /i-dependent energy eigenvalues of Hg (/i) , Es^k (m): 



Ek < minEs.kil^) , k 

M>0 



0,1,2,. 



For a pure Coulomb potential y(r) — —a/r, the energy 
eigenvalues i?s.„(/i) of iJg (/i) read, with rt = n^+i+l, 



2 II 



1-^ 



Here, minimizing £^s.n(M) with respect to /i entails [17] 



uiinEs.nil^) = TO Wl - 
M>o ■ V 

This (exact) upper bound [ 17] to the energy eigenvalues 
of the so-called "spinless relativistic Coulomb problem" 
holds for all those values of the Coulomb coupling a for 
which the Hamiltonian H with a Coulomb potential can 
be regarded as a reasonable operator and arbitrary levels 
of excitation, and for any value of the principal quantum 
number n it definitely improves the Schrodinger bound: 

min£'s,n(M) < Es,n for a ^ . 
Clearly, fixing = m recovers the Schrodinger bounds. 
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3.4 Exact semianalytical upper and lower 
bounds from the "envelope technique" 

Rigorous semianalytical expressions for both upper and 
lower bounds to the eigenvalues Ene of the Hamiltonian 
H are found by a geometrical operator comparison in an 
approach called "envelope theory." The envelope theory 
constructs bounds on Eni by comparing the spectrum of 
H with the one of a conveniently formulated "tangential 
Hamiltonian" H involving some "basis potential" h{r), 



H = \/mM-p^ + c /i(r) + const. , c > , 

for which sufficient spectral information (i. e., either the 
exact eigenvalues or suitable bounds on these) is known. 
Let V{r) be a smooth transformation V = g{h) of h{r), 
with definite convexity of g{h). After optimization with 
respect to the point of contact of y (r) and the tangential 
potential, this technique produces bounds on Ene'- lower 
bounds for g{h) convex (g" > 0) , and upper bounds for 
g{h) concave {g" < 0). Suppressing for the moment the 
quantum numbers n£, all these bounds on E may be cast 
into a common generic form [20-25] with the individual 
bounds discriminated by a dimensionless parameter, P: 



E' 



mm 

r>0 



(12) 



Here, that cryptic sign of approximate equality indicates 
that for any definite convexity of g{h) all expressions on 
the right-hand side represent a lower bound for a convex 
g{h) and an upper bound for a concave g{h). The value 
of the parameter P used in Eq. (I12> is determined by the 
algebraic structure of the interaction potential V{r), and 
by its convexity with respect to the basis potential, h{r): 

• The spinless relativistic Coulomb problem posed 
by V{r) = — a/r is well-defined if its coupling a 
is constrained to a < 2/-K [26] . The bottom of the 
corresponding spectrum of H (or, its ground-state 
energy eigenvalue, Eq) is bounded from below by 



EQ>m\ \ 



p2 ' 



with the lower-bound parameter P given either by 
P = 2/7r for a fulfilling < a < 2/7r [26], or by 



P = P{a) 



1 + -4( 



for < a < ^ , which obviously covers the range 

P{\) = ^ < P{^) < P{0) = 1 , 

as derived by weakening [21-23,25] an improved 
lower bound to Eq valid only for < a < ^ [27] . 



If V{r) is a convex transform V ~ g{h), g" > 0, 
ofthe Coulomb potential /i(r) = —1/r, the above 
envelope approximation generates a lower bound 
[21-23,25] on the ground-state eigenvalue Eq (on 
the entire spectrum) of the Hamiltonian H for any 
choice of the Coulomb lower-bound parameter P. 

Clearly, the quoted upper bounds on the Coulomb 
coupling a apply also to any "effective" Coulomb 
coupling in H. Consequently, they translate into a 
constraint on all coupling constants introduced by 
the interaction potential V{r) under investigation. 
(An example for these restrictions enforced by the 
Coulomb menace will be given in Subsect. r3.7.2l ) 

For V{r) a concave transform V ~ g{h), g" < 0, 
of the harmonic-oscillator potential h{r) — r^, a 
straightforward application of the above envelope 
approximation yields upperhonnds [21-23,25] to 
all the eigenvalues of the Hamiltonian H; the 
parameter P for a given energy level identified by 
quantum numbers n£ is, in this case, related to the 
explicitly algebraically known eigenvalues of 
the (nonrelativistic) Schrodinger operator p^+r^: 



P = Pne{2)^^£ne = 2n + £- 



FoiV{r) a concave transform y = g{h), g" < 0, 
of the linear potential h{r) = r, the application of 
a "generalized" envelope approximation provides 
M/jper bounds [23,25] to all eigenvalues Eni of H 
if the parameters P which characterize the energy 
levels are given, in terms of the eigenvalues of 
the nonrelativistic Schrodinger operator p^+r, by 



P = P„,(l) = 2(i£„, 



,3/2 



(13) 



the parameter values corresponding to the 

lowest-lying energy levels can be found in Table[2 
(for more details see, for instance, Refs. [20-23]). 

If the potential V{r) is the sum of several distinct terms, 

i 

where every component problem defined by the operator 



supports, for a sufficiently large Ci, a discrete eigenvalue 
i?i.o at the bottom of its spectrum and information about 
the lowest energy eigenvalue, E'i.o, is available, all these 
pieces of information can be combined to a lower bound 
to Eq [24] ; for sums of pure power-law terms sgn(g) , 



VpL(r) = ^ a{q) sgn(g) r« 



(14) 
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Table 1 : Numerical values of the parameter Pne{l) used 
in the linear-potential-based lower envelope bounds and 
defined in Eq. il3\ for the lowest-lying energy levels n£. 



n 


£ 


P„f(l) 


n 


i 




1 





1.37608 


1 


3 


4.36923 


z 


A 
U 


1 8 1 1 


z 




0. i jZVo 


J 


A 
U 




J 


J 


7 Ql '^flA 


A 
H- 


A 
U 


f. QA« 1 A 


A 


J 


v. /UZjO 




A 
U 


o.OioZJ 




J 




1 


1 


2.37192 


1 


4 


5.36863 


2 


1 


4.15501 


2 


4 


7.12732 


3 


1 


5.95300 


3 


4 


8.90148 


4 


1 


7.75701 


4 


4 


10.68521 


5 


1 


9.56408 


5 


4 


12.47532 


1 


2 


3.37018 


1 


5 


6.36822 


2 


2 


5.14135 


2 


5 


8.12324 


3 


2 


6.92911 


3 


5 


9.89276 


4 


2 


8.72515 


4 


5 


11.67183 


5 


2 


10.52596 


5 


5 


13.45756 



where the coefficients a{q) of the pure power-law terms, 
sgn(g) r'^ , in the potential are positive, that is, a{q) > 0, 
and do not vanish all, this yields the "sum lower bound" 



Eq > min 

r>0 



^ + J2ciiq) Sgn(q) (P(g) r)« 



provided that some set of lower-bound parameters P,{q) 
can be derived such that, whenever V{r) consists of just 
one single component, the above inequality yields either 
the corresponding exact ground-state energy eigenvalue 
or, at least, a rigorous lower bound to this latter quantity: 

• For Coulomb components, that is, hi (r) = — 1 /r, 
P(— 1) is the Coulomb lower-bound parameter P. 

• For linear components, that is, hi{r) ~ r, P_{l)is 
derived from the lowest eigenvalue Sq of i/p^+r, 



£(1) 



1.2457. 



It is straightforward to (try to) generalize these envelope 
techniques from the simpler one-body case summarized 
in this review to systems composed of arbitrary numbers 
of relativistically moving interacting particles described 
by a semirelativistic spinless Salpeter equation [28-30]. 
At least for the particular case of all harmonic-oscillator 



bounds presented in Subsect. l33l and the envelope upper 
bounds can be shown to be equivalent to each other [23] . 

3.5 Rayleigh-Ritz (variational) technique 

An immediate consequence of the minimum-maximum 
principle is the "Rayleigh-Ritz (variational) technique:" 

• Introduce the restriction H of some operator H to 
a subspace Dd by orthogonal projection P to Dd'. 



H = H 



:= PHP 



• Identify all d eigenvalues Ek, k = 0,1, . . . , rf— 1, 
of the restricted operator H as the solutions of the 
eigenvalue equation of H for the eigenstates \xk)- 



H\xk)^Ek\xk) , A: = 0,1, 



1 



• Let these eigenvalues Ek be ordered, according to 

Eq < El < ■ ■ ■ < Ed^i . 

Then every (discrete) eigenvalue Ek of H — if counting 
the multiplicity of degenerate levels — is bounded from 
above by the eigenvalue Ek of the restricted operator H: 



Ek<Ek for all fc = 0.1, 



1 



If that d-dimensional subspace Dd is spanned by any set 
of d (of course, linearly independent) basis vectors \ipk), 
k = 0,1, . . . , 1, the eigenvalues Ek can immediately 
be determined, by the diagonalization of the dxd matrix 



potentials V{r) 



with c > the generalized upper 



{^,\H\^,)j, i,j^O,l,...,d-l, 

that is, as the d roots of the characteristic equation of H, 

det((^,|#|V,)-^(V'^lV',>) = 0, 
«,J = 0, 1, . . . ,d - 1 . 

To establish this, expand any eigenvector \xk) of H over 
the basis i — 0,1, ... , d—1} of the subspace Dd. 

3.6 Variational upper bounds 

The quality achieved by the variational solution of some 
eigenvalue problem depends decisively on the definition 
of the trial subspace Dd employed by the Rayleigh-Ritz 
technique briefly sketched in Subsect. 1331 enlarging Dd 
to higher dimensions d or choosing a more sophisticated 
basis {It/ji), i = 0,1, . . . , d—1} which spans Dd will, in 
general, increase the accuracy of the obtained solutions. 
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For spherically symmetric (central) potentials V{r), 
that is, for all potentials which depend only on the radial 
coordinate r = \x\,a convenient and thus rather popular 
choice for the basis vectors {iV'i), i = 0, 1, . . . , d~l} is 
that one the configuration-space representation of which 
involves the complete orthogonal system of generalized 
Laguerre polynomials [12,13,31,32] — cf. AppendixlAl 



In the one-dimensional case [33] realized in the notation 
of AppendixElif all quantum numbers k = £ = m = 0, 
the Laguerre basis collapses to just a single basis vector: 



exp(— /irj 



/I > 



With a trial state lip) represented by this exponential and 
the trivial (nevertheless fundamental) general inequality 



mom ^ 



which holds for any self-adjoint, but otherwise arbitrary, 
operator© (O^^ — O), Rayleigh's principle entails, after 
optimization with respect to the variational parameter fi, 
for a Coulomb potential V{r) = —a jr the upper bound 



Eq < m \/\ — cP- ; 

this is identical to the "generalized" upper energy bound 
on the ground-state or n = 1 eigenvalue of the Coulomb 
operator H found by different reasoning in Subsect. l33l 



3.7 Application to illustrative interactions 

Let us appreciate the above bounds' beauty at examples. 



3.7.1 Trivial "testing ground:" Coulomb potential 

Our first example clearly must be the Coulomb potential 

y(r)=--, a>0; 
r 

this potential arises from the exchange of some massless 
boson between the interacting objects. Therefore it is of 
particular interest in many areas of physics. Its effective 
interaction strength is given by a coupling a, identical to 
the fine structure constant in electrodynamics. We study 

• the somewhat naive nonrelativistic (Schrodinger) 
upper bound given in Subsect. r3.2.1l equivalent to 
a tangent line to the relativistic kinetic operator T, 

• the upper bound of Subsect. r3.2.2l constructed by 
considering just the square of the Hamiltonian i/, 

• the semianalytical upper bound of Subsect. l33l as 
derived by generalizing the idea of Subsect. r3.2.1l 

• all three envelope bounds of Subsect. l3r4l namely, 

- the harmonic-oscillator-based upper bound, 

- the upper bound involving a linear potential, 

- the lower bound obtained by "loosening" an 
absolute lower bound on the spectrum of the 
"semirelativistic Coulomb operator" i7, and 

• the "Rayleigh-Ritz" upper bound of Subsect. ll!6l 
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For the Coulomb potential under study, the optimization 
required by the envelope bounds (I12> may be performed 
analytically, yielding a result of precisely the same form 
as the generalized upper bounds derived in Subsect. l33l 
or as the squared upper bounds proved in Subsect. r3.2.2l 



Eo{P) 



1 



P2 ' 



(15) 



where for the ground state characterized by the quantum 
numbers n — 1,£ — Q the (single) parameter P is given. 



• for the "Coulomb lower bound" (Subsect. l3.4t . by 



c 



Pia) 



• for the generaUzed upper bound (Subsect. l3.3> . by 

P = PG=n = l, 



• for the "linear upper bound" (Subsect. l3.4> . as can 
be simply read off from the first row in Table[T] by 

P = Pl=Ao(1) = 1.37608, 



• for the "squared upper bound" (Subsect. l3.2.2t . in 
accordance with the solution of Eq. {HJ for L, by 



P = Po 



V2N 



1 + VT+Jo? 

V2 



• and, in the case of the "harmonic-oscillator upper 
bound" (Subsect. l3.4> . from the Pni (2) results, by 



Ao(2) 



It is a very trivial observation that, for fixed values of the 
Coulomb coupling, the ground-state energy bounds il5\ 
are (monotone) increasing with increasing parameter P: 



dP 



> . 



Thus it is straightforward to convince oneself that all the 
Coulomb (Ec), generalized (Eg), nonrelativistic (-En), 
linear (E'l), squared (Eq) and harmonic-oscillator (£'h) 
bounds on the ground-state energy eigenvalue Eq of the 
semirelativistic Coulomb Hamiltonian H have to satisfy 

Ec ^ Eq < Eq < i?N < -Bl < Eq < E^l 

for a<ao= f (s - 2 \/2) , 
Eq < Eq < Eq < £'n < -El < -Eh < Eq 
fora > uq = W| fs- 2\/2 



taking into account the crossing of the upper bounds Eh 
andEq at = | (3-2 \/2), i. e., Eq{ao) = EH(ao). 



For Coulomb-like interactions the only dimensional 
quantity among the parameters of this theory is the mass 
m of the interacting particles. Consequently, in this case 
all energy eigenvalues are proportional to m: the energy 
scale is set by m. The ratio E^ /m is a universal function 
of the coupling a; w. 1. o. g. it thus suffices to fix m = 1. 

Figure^compares for the ground state (rii — I — 0) 
of the spinless relativistic Coulomb problem the various 
bounds to the lowest energy eigenvalue, Eq , listed at the 
beginning of this subsection. Inspecting Fig.^ we note: 

• the squared, harmonic-oscillator, and linear upper 
bounds are numerically comparable to each other; 

• likewise the nonrelativistic and generalized upper 
bounds are close to each other for all couplings a; 

• using a Laguerre trial space of dimension d = 25, 
the Rayleigh-Ritz variational upper bound can be 
expected to come akeady pretty close to the exact 
eigenvalue Eq — which, in turn, clearly indicates 
that it is highly desirable to find improvements for 
the lower bounds, in particular for large couplings 
a (this stimulated, e. g., the analysis of Ref. [34]). 



0.95 



E 0.85 




0.75 



Figure 1 : Both upper (full lines) and lower (dashed line) 
bounds on the ground-state energy eigenvalue (E) of the 
semirelativistic Hamiltonian H with Coulomb potential 
V{r) = —a/r as a function of the Coulomb coupling, a, 
for the squared (Q), harmonic-oscillator (H), linear (L), 
nonrelativistic (N), generalized (G), variational (V) and 
Coulomb (C) [using the "optimized" P{a)] approaches. 
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3.7.2 Coulomb-plus-linear (or "funnel") potential 

Within the field of elementary particle physics, quantum 
chromodynamics (QCD) is generally accepted to be that 
relativistic quantum field theory that describes all strong 
interactions between quarks and gluons by assigning the 
so-called "colour" degrees of freedom to these particles. 
In the instantaneous approximation inherent to all of the 
QCD-inspired quark potential models developed for the 
purely phenomenological description of experimentally 
observed hadrons, as bound states of constituent quarks, 
the strong forces are assumed to derive from an effective 
potential generating the bound states (this description of 
hadrons within the framework of quark potential models 
involving either nonrelativistic or relativistic kinematics 
isreviewed, for instance, in Refs. [3,35].) The prototype 
of all "realistic," that is, phenomenologically acceptable 
(static) interquark potentials V{r) consists of the sum of 

• a Coulomb contribution generated by a one-gluon 
exchange between quark bound-state constituents 
(dominating the potential at short distances r) and 

• a linear term including all nonperturbative effects 
(that dominates the potential at large distances r). 

The resulting interaction potential V{r) is characterized 
by a "funnel-type" Coulomb-plus-linear form; therefore 
it is called the Coulomb-plus-linear, or funnel, potential. 
Upon factorizing off a constant v, which spans the range 
< u < 1 in order to parametrize an overall interaction 
strength, we (prefer to) analyze this potential in the form 

V{r) = -— + C2r = v(-- + br) . (16) 
r \ r J 

Clearly, given the overall coupling strength u, the actual 

shape of this potential is fixed by the ratio of the positive 

parameters a > and 6 > 0; the coupling constants that 

enter, on the one hand, in the general expression MAI for 

sums of pure power-law terms and, on the other hand, in 

our funnel potential M6\ must be identified according to 

a(— 1) = ci = a w > , 

a(l) = C2 = hv>Q. 

In view of the lack of fully analytical bounds we explore 

• the three "basic" envelope bounds of Subsect. 1341 
distinguished by the adopted basis potential, viz., 

- the upper bound from a harmonic oscillator, 

- the upper bound involving a linear potential, 

- the lower bound due to a Coulomb potential, 

• the envelope sum lower bound, derived in the sum 
approximation recalled by Subsect. lT?! as well as 

• the "Rayleigh-Ritz" upper bound of Subsect. ll!6l 



For definiteness, let us fix the potential parameters a and 
6 to a = 0.2, h = 0.5. As done in the Coulomb-potential 
example (in Subsect. r3.7.1> in order to take advantage of 
upper and lower bounds, we investigate the ground-state 
energy E'q. The basic envelope bounds are computed by 
application of Eq. (I12> . for the appropriate parameter P: 

• for the "harmonic-oscillator upper bound" we use 

P = Ph = Pio(2) = I ; 

• for the "linear upper bound" we find from Table[2 

P = Pl=Ao(1) = 1.37608; 

• for the "Coulomb lower bound," in order to derive 
the maximum value P consistent with < w < 1 , 
we are forced to evaluate that "Coulomb coupling 
constant constraint" mentioned in Subsect.|^3 in 
its form [22,23] fixed by our funnel potential J16L 

for the maximum values of ci and C2 , which gives 
P = Pc = 0.728112397 for m = 1 . 

The "sum lower bound" is extracted from the expression 
given explicitly in Subsect. l3r4l for power-law potentials 
by insertion of the lower-bound parameters P_{q = ±1): 

• the Coulomb lower-bound parameter P(a) leads, 
for the relevant maximum coupling a = a = 0.2, 
in the Coulomb term of the sum approximation to 

P(-l) = P(a) = P(a) = 0.9789063 ; 

• the lower-bound parameter required for any linear 
part of sum potentials is copied from Subsect. l3r4l 

P(l) 1.2457. 

As before, the Rayleigh-Ritz or variational upper bound 
is found in a trial space of dimension d = 25 spanned by 
the generalized Laguerre basis (summarized in App.lAl. 

FigurelJldepicts the bounds to the lowest eigenvalue 
Eq of H as function of the overall coupling strength v in 
the funnel potential il6i . Remarkably, variational upper 
and sum lower bounds now restrict Eq to a narrow band. 

4 Approximate Solutions: Quality 

Having determined — for instance, by application of the 
Rayleigh-Ritz technique sketched in Subsect. l331 — for 
some fc = 0, 1, 2, . . . the state \xk) G Dd corresponding 
to any upper bound Ek on the exact eigenvalue Ek of H, 
one question immediately arises: how closely resembles 
the approximate solution \xk) the exact eigenstate ? 
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Figure 2: Three upper (fuJJ Jines) and two lower (dashed 
lines) bounds on the ground-state energy eigenvalue (E) 
of the semirelativistic Hamiltonian H with the so-called 
funnel potential V{r) = v {—a/r + br), where a = 0.2, 
b — 0.5, m — 1. These include: the harmonic-oscillator 
(H), linear (L) and variational (V) upper bounds and the 
sum-approximation (S) and Coulomb (C) lower bounds. 



Standard criteria, such as the (relative) distance between 
Ek and true Ek , or the overlap of approximate and exact 
eigenstates, require the knowledge of the exact solution. 
In contrast to this, the virial theorem (Sect.|2j represents 
an indicator for the accuracy of approximate eigenstates 
that merely uses information provided by the variational 
approach: Since all eigenstates of H satisfy any relation 
of the form (|5}, a significant imbalance in Eq. Q reveals 
that this approximation is far from optimum [23, 36, 37] . 
Of course, because of the involvement Q of the dilation 
generator (|6j in the derivation of Eq. (0, any variational 
solution found by minimization of expectation values of 
H with respect to the scale transformations, or dilations, 
Q will necessarily satisfy our master virial theorem 0. 



5 Summary, Concluding Remarks 

The various efficient approaches presented here allow to 
analyze the semirelativistic Hamiltonians of the spinless 
Salpeter equation analytically; this is crucial for general 
considerations that aim to answer questions of principle, 
like operator boundedness. For numerical methods, see, 
for instance, Refs. [10,38,39] and the references therein. 



A The Generalized Laguerre Basis 

Assume every basis function of L2 (i?"^ ) to factorize into 
a function of the radial variable and the angular term. Its 
configuration-space representation has the general form 

^k,ir,i{x) = '^kA^)yirn{^r) , r = \x\ ] 

the spherical harmonics (fi) for angular momentum 
I and projection m depend on the solid angle Q. = {9,(j)) 
and satisfy a well-known orthonormalization condition: 



The most popular choice [12,13,31,32] for the basis 
states which span the Hilbert space L2{R^) of [with the 
weight w(a;) — x^] square-integrable functions f{x) on 
the positive real line i?+ — which is the Hilbert space of 
radial trial functions ^k.e (r) — involves the generalized 
Laguerre polynomials l'"^\x), for parameter 7 [40,41]: 



these generalized Laguerre polynomials for parameter 7 
are orthogonal polynomials, defined by the power series 



k + J 
k ~ t 



X' 



0,1, 



t=0 

and orthonormalized with weight function x'^ exp(— x): 

00 

J dxa;^exp(-x)4^'(x)4T^(x) 


r(7 + fc + i) 



ki 



^kk' 



k,k' = 0,1,, 



The basis states defined by the generalized-Laguerre 
choice for the radial basis functions ^k.i (r) involve two 
parameters, both of which may be subsequently adopted 
for variational purposes: ^ (with the dimension of mass) 
and f3 (dimensionless); requirements of normalizability 
of our basis states constrain the parameters to the ranges 

< ^ < 00 , -1<2/3<(X). 

Therein, the orthonormality of the generalized Laguerre 
polynomials, inherent to their definition, is equivalent to 
the orthonormality of the radial basis functions ^k,e{r)' 



drr^^kAr)'^k'Ar)=Skk' , k, k' = 0,1, . . . ; 



this condition fixes the normalization constant Nj.'^^^^ to 



k,e 



T{2e + 2(3 + k + l) 
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Fortunately the assumed factorization of every basis 
function persists in its momentum-space representation: 



Analytical statements about Hamiltonians that involve a 
kinetic-energy operator nonlocal in configuration space, 
such as a relativistic square root Q, are facilitated by an 
explicit knowledge of the momentum-space basis states. 
One of the great advantages of the generalized-Laguerre 
basis is the availability of its analytic Fourier transform. 

For all factorizations into radial and angular parts, as 
consequence of the Fourier transformation acting on the 
Hilbert space L2 (R^) of the square-integrable functions 
on the three-dimensional space R^, the radial parts of all 
basis functions that represent the chosen basis vectors in 
configuration space and momentum space, respectively, 
are related by so-called Fourier-Bessel transformations: 



^kAp) = {-V y - J <irr^ je{pr)<!>kAr) , 


forallfc = 0, 1,..., ^ = 0, 1,... ; 

the angular-integration remnants jn{z) {n ~ 0, ±1, . . .) 
label the spherical Bessel functions of the first kind [40]. 
For the generalized-Laguerre basis under consideration, 
these radial basis functions become in momentum space 



k,e 



{—iYp^ 



2^+1/2 r(£ + 1) 

(-1)* / k + 2£ + 2f3 



t=o 



k - t 



2 ' 2 ' 2>2 + ^2y' 

with the hypergeometric series v] w; z), defined, in 
terms of the gamma function F, by the power series [40] 

F{u, v; w; z) 

T{w) ^ r{u + n)T{v + n) z"" 
~ r(u)r(v) ^ T(w + n) ^ ' 

and the simplifying abbreviation at. f;/3 = 2 £+ /3+t + 2. 

Clearly, the momentum-space radial basis functions 
^k,e{p) have to satisfy the orthonormalization condition 



kk' 



k,k' = 0,1,. 
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